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Based on very accurate measurements performed on a superconducting microwave resonator 
shaped like a desymmetrized three-dimensional (3D) Sinai billiard, we investigate for the first time 
spectral properties of the vectorial Helmholtz, i.e. non-quantum wave equation for a classically 
totally chaotic and theoretically precisely studied system. We are thereby able to generalize some 
aspects of quantum chaos and present some results which are consequences of the polarization 
features of the electromagnetic waves. 
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For nearly 20 years, billiard systems have provided a 
very effective tool for the investigation of semiclassical 
quantization of conservative chaotic systems (l[|^|. This 
is due to the fact that even two-dimensional (2D) bil¬ 
liards (as opposed to billiards of higher dimensionality) 
are able to model a wide range of fully ergodic systems 
in Gutzwiller’s sense of ’’hard chaos” |j|. As a matter of 
fact, properties of the wave dynamical spectra of such 
low-dimensional but classically non-integrable systems 
are fully described by the Gaussian Orthogonal Ensem¬ 
ble (GOE) of Random Matrix Theory (RMT) |||| if the 
underlying motion is invariant under time-reversal. On 
the other hand, classically regular, i.e. integrable systems 
lead to totally uncorrelated spectra. 

Up to now investigations on chaotic 3D-billiards were 
performed in experiments with electromagnetic j6|7]] and 
acoustic waves, whereas the hardly feasible numeri¬ 
cal modelling was restricted to very special geometries of 
high symmetry for the pure Schrodinger problem |l(j|j . 
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FIG. 1. Geometry of the desymmetrized 3D-Sinai billiard 
(boldface line) which constitutes one sixth of the dashed cube. 
Eight of those cubes form the full system. As indicated in the 
figure, one antenna was located in the center of each plane 
surface of the cavity. 


The goal of the present paper is to provide for the 
first time a detailed analysis of a fully chaotic three¬ 


dimensional (3D) electromagnetic billiard with a classi¬ 
cally well-known and theoretically precisely studied ge¬ 
ometry: the 3D-Sinai billiard resp its desymmetrized 
version given by 1/48 of a cube with a sphere in its 
center, see Fig.l. The system has to be described by 
the time-independent, fully vectorial Helmholtz equation 
with electromagnetic boundary conditions. The same 
geometry was recently investigated numerically in the 
quantum regime described by the time-independent 
Schrodinger equation and experimentally with acoustic 
waves [^|. Our analysis therefore allows a very distinct 
comparison of totally different wave dynamical phenom¬ 
ena in a system with exactly the same classical analogue. 
Our results should agree with those of |u| and Q if the 
conjecture holds that RMT is adequate to describe spec¬ 
tra of arbitrary wave phenomena. 

As several 2D- and 3D-billiards before iHMl , the 
electromagnetic resonator was made of Niobium which 
becomes superconducting below 9.2 K. This feature 
tremendously increases the resolution of the measured 
spectra due to quality factors of up to 10 7 compared to 
10 3 in normal conducting resonators and signal-to-noise 
ratios of up to 70 dB. The resonator was cooled down to a 
temperature of 4.2 K and finally evacuated to a pressure 
of 0.2 mbar inside a bath cryostat. We were able to excite 
the cavity via four different antennas up to a frequency 
limit of 20 GHz given by our HP8510B vector network 
analyzer which was used for the generation and detec¬ 
tion of the microwave signal. Either the reflected signal 
from one antenna, or the transmitted signal between two 
different antennas was measured. Figure 2 shows a typ¬ 
ical transmission spectrum between 6.50 and 6.75 GHz. 
Identifying accurately per hand the positions of signifi¬ 
cant peaks of every single spectrum and performing in 
the next step a detailed comparison of all the possible 
spectra, a total set of 1881 experimental eigenmodes was 
obtained which form the base of the following investiga¬ 
tions. 
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To check for the mechanical accuracy of the cavity and 
to estimate the effects due to the coupling of the iso¬ 
lated resonator to the external world (coupling holes and 
slightly penetrating coupling wires of the rf-cables), we 
also simulated the spectrum using the electromagnetic 
CAD-program MAFIA (solution of Maxwell’s equation 
by the Finite Integration Algorithm [Q). The resonator 
was modelled using a mesh consisting of 10 6 points, and 
a full 3D broadband time domain calculation with a to¬ 
tal number of 150,000 time steps was performed. A task 
which is feasible on a standard workstation with one week 
of CPU time. The eigenfrequencies were finally deter¬ 
mined by a Fast Fourier Transform yielding a resolution 
of about 1.2 MHz over a wide frequency range. Limited 
by this resolution a comparison with the experimental 
spectrum (with a resolution of 10 kHz) allows a one-to- 
one correspondence between simulation and experiment 
within a region up to 7.5 GHz including 94 resonances. 

For comparison, the numerically simulated eigenfre¬ 
quencies are given as dashed vertical lines in Fig.2. As 
a matter of fact the experimental eigenfrequencies are 
systematically up-shifted due to the contraction of the 
resonator during the cool-down. From this shift (linearly 
increasing with frequency) we extracted the relative con¬ 
traction coefficient ( Ax/x) exp ~ 2.0 • 10~ 3 which is close 
to the value for Niobium as given in the literature [fb5| , 
(A x/x)ut ~ 1.4 • 10 -3 . Beside this effect, every experi¬ 
mental resonance is individually shifted due to the cou¬ 
pling holes and antennas and the mechanical fabrication 
tolerances. However, these individual shifts are on the 
order of some 5 • 10 “ 4 and thus much smaller than the 
effect of the cool-down. 

To prepare the experimental spectrum of extracted 
eigenfrequencies for the statistical analysis, we first un¬ 
folded it, i.e. rescaled the frequency axis to a mean level 
spacing of unity. According to this, in applying the gen¬ 


eralized electromagnetic Weyl formula 16 1 
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to describe the smooth part of the integrated level den¬ 
sity (spectral staircase), the spectrum is again checked for 
completeness. Here, N counts the total number of reso¬ 
nances up to a certain frequency /, V denotes the volume 
of the geometry and Co is the speed of light. In the linear 
term, R labels the mean radius of curvature over the sur¬ 
face cr and is the dihedral angle along the edges a. Note 
that the leading cubic term of this expression is twice the 
corresponding term in the scalar, i.e. quantum problem, 
and that there is no quadratic term . Both features are 
due to the vectorial character of the electromagnetic field 


inside the cavity and thus a consequence of two trans¬ 
verse polarizations of the electromagnetic wave. A com- 
parision between Eq. ®) and the integrated experimental 
level density yields an estimate for the number of uncer¬ 
tain resonances AN ss ( V exp /V t heo - 1)A r totai ~ 4. Here, 
Aftotai denotes the complete set of 1881 resonances. Af¬ 
ter the extraction of this smooth part the expected non- 
systematic fluctuations around zero could be observed, 
they are the carrier of all accessible information about 
the classical dynamics inside the billiard. 
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FIG. 2. Excerpt of a typical transmission spectrum. The 
signal displays the ratio of output power to input power on a 
logarithmic scale. The numerically simulated eigenfrequencies 
are given by dashed lines. The shift between experiment and 
simulation is mainly due to the contraction of the resonator 
during the cool-down. 

Statistical measures developed to display the correla¬ 
tions embedded in these fluctuations were originally in¬ 
troduced in the 1950th to describe spectra of nuclear sys¬ 
tems with many degrees of freedom. Q. An adequate 
and widely tested statistical measure to examine short- 
range correlations up to a length of one mean level spac¬ 
ing is given by the nearest neighbour spacing distribu¬ 
tion (NND), displaying the probability P(s) for a certain 
spacing s of two adjacent resonances in the spectrum. 
Figure 3 (l.h.s.) displays P(s) of the experimental spec¬ 
trum as well as the two limiting curves, ’’Poisson” for a 
totally uncorrelated spectrum corresponding to regular 
motion on the classical side, and ”GOE” for the totally 
chaotic case. Obviously, the experimental histogram is 
very close to the GOE distribution. In the next step 
we analyzed the spectrum on a larger scale in order to 
investigate long-range correlations. For this purpose we 
calculated E 2 (L), which expresses the variance of a num¬ 
ber of resonances inside an interval of length L on the 
unfolded scale, as well as the related Dyson-Mehta statis¬ 
tics, A 3 (L), also sensitive up to L, i.e. several mean level 
spacings. The result for both properties is given as well in 
Fig.3 (l.h.s.). Here, two observations can be made: First, 
the experimental curves rapidly deviate from the GOE 
prediction and lie in between the regular and the chaotic 
case, and second, above a certain value L max , which is 
different for both statistics (L^ al « 40, L^ ax « 150), 
the experimental curves run into saturation. This last 


2 





























feature is exactly what is expected from theory [|T8j| , dis¬ 
playing the fact that for increasing L the given statistics 
are more and more sensitive to specific, i.e. non-universal 
features of the system. In this sense L max separates the 
universal scale which is dominated by the general dy¬ 
namical behaviour of the classical analogue - described 
by the collective behaviour of all long periodic orbits of 
the system - from the non-universal scale which is gen¬ 
erated by individual system-specific attributes - in the 
present case given by the shortest periodic orbits of the 
billiard. Note that the difference between L max of E 2 
and the value obtained for A 3 is t^ 3 OJ « 4 • L^ ax [^9| , 
which is well reproduced by the given spectrum. The¬ 
oretically the length of the shortest periodic orbits l m in 
depends on L%* ax according to 0 
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where f max denotes the upper frequency limit of 20 GHz. 
Using the length of the first observable periodic orbit, 
lmin= 0-34 m (see below), one estimates L^ ax ss 125, 
which is close to the value of Fig.3. 
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FIG. 3. Short-range and long-range statistical measures 
before and after the extraction of 3D bbos, see text. 


To characterize all statistical measures more quantita¬ 
tively and in order to check how well the GOE prediction 
describes the results, we also analyzed the experimental 
curves according to a model of Berry and Robnik J2(]] 
within the universal region. The basic intention of this 
model is to provide a continuous interpolation between 
the pure Poissonian and GOE using a mixing-parameter 
q which corresponds to the relative phase space volume 
that is covered by chaotic trajectories. In the present 
case we obtained for the three statistical measures on 
the l.h.s. of Fig.3 qNND = 0.95 ± 0.01, 2 = 0.89 ± 0.02 

and < 7 a 3 = 0.91 ± 0.02, indicating a small but systematic 
deviation from unity, the limit for pure chaotic systems. 
Thus, although the system is fully ergodic on the classical 


side without any stable islands in phase space, the wave 
dynamical side pretends their existence. To understand 
this phenomenon, we analyzed the classical analogue in 
more detail. Therefore, we continued our investigation 
on a more specific scale which provides a bridge between 
the classical chaotic features and their impact on the elec¬ 
tromagnetic spectrum. This scale is given by the length 
spectrum of the billiard, which is practically obtained 
through the Fourier transform of the spectral level den¬ 
sity or its fluctuating part, respectively, 

P fmax p_ 

P fluc (l ) = / P fluc (f) ■ exp(*— If) df . (3) 

Jfmin C 0 

Here l denotes the physical length scale. Since the prop¬ 
erty pf luc (l ) maps the long periods in frequency space 
onto short scales in length space, the resulting spectrum 
shows peaks near the classical periodic orbits of the bil¬ 
liard. Figure 4 exhibits the lower part of this length spec¬ 
trum up to Z=1.5 nr. Here a rich structure of peaks can 
be observed above a minimum length Z m m=0.34 m. This 
first peak belongs to the shortest 3D bouncing ball orbit 
(3D bbo) of the billiard, propagating along one edge of 
the desymmetrized cube without striking the sphere, see 
sketch in the inset. 



Orbit length l (m) 

FIG. 4. Experimental length spectrum of the billiard (full 
line) and the semiclassical reconstruction using only 3D bbos 
(dashed line). The pictures in the inset show the first unstable 
and the first 3D bouncing ball orbit, respectively. 

Exactly this shortest length was used to estimate the 
saturation parameter L max of A 3 (see above). In fact, 
this first 3D bbo is the shortest non-isolated, neutrally 
stable periodic orbit of the system. It is highly remark¬ 
able, that this is not the first periodic orbit of the bil¬ 
liard. The shortest unstable periodic orbit runs along the 
shortest edge of the desymmetrized billiard which is inter¬ 
sected by the sphere, it possesses a length l^° n =0.17 m, 
see sketch in the inset. The corresponding peak in the 
Fourier spectrum is drastically suppressed due to the 
desymmetrization of the system itself |To],^l]] . As a mat¬ 
ter of fact, the obtained length spectrum is totally domi¬ 
nated by bbos of all possible dimensions, not only in the 
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quantum case |T(J but also in the electromagnetic coun¬ 
terpart. To demonstrate this, we considered the contri¬ 
bution of the leading 3D bbos to the given length spec¬ 
trum. Therefore, we used a lattice vector description 
of Berry [Q to label and generate all 3D bbos up to 
a certain length and determined their contribution on 
p{l U oV) = dN^ c /df through @ 

nSS(X) = ^ ( •£ (x 2 - „’) - Ar» + V) , 

l bbo \ 0<n<X ^ Z ' 

( 4 ) 

with X = Ibbof/co . Here, the length of a given 3D 
bbo, lbbo-, was deduced directly from the lattice vector 
and Sbbo , the perpendicular area on which this orbit ex¬ 
ists, was fixed in a Monte-Carlo simulation. In the given 
range up to Z=1.5 m we obtained 55 3D bbos of differ¬ 
ent degeneracies and with positive Sbbo-, their superim¬ 
posed semiclassical reconstruction due to Eq. 0) is given 
by the dashed line in Fig.4. It is highly remarkable that 
nearly the full structure of the given length spectrum can 
be reproduced using only 3D bbos, whereas the influ¬ 
ence of the enormous number of unstable periodic orbits 
(approximately 36000 up to 1= 1.5 m HD is hidden in 
the background. Discrepancies between the experimental 
length spectrum and the reconstruction can be predomi¬ 
nantly found at the locations of the 3D bbos themselves 
and arise because of the existence of sub-dimensional and 
tangential bbo manifolds (n]] . 

Finally, to demonstrate the influence of the considered 
3D bbos on the statistical measures NND, E 2 and A 3 , we 
repeated our statistical analysis using a modified unfold¬ 
ing procedure in which the standard Weylian, Eq. ( 0 ) , is 
extracted together with the contribution of the 3D bbos, 
Eq.(||). For this we used all 3D bbos up to a length 
Z=3.0 m. The result is given in Fig.3 (r.h.s.) displaying, 
now also for E 2 and A 3 , nearly perfect agreement with 
the GOE prediction in the universal regime up to L max . 
The corresponding, i.e. corrected mixing parameters, 
were fixed to be qNND = 0.96 ± 0.01, 2 = 0.99 ± 0.01 

and (?a 3 = 0.99 ± 0.01. 

In summary, a set of 1881 highly resolved eigenmodes 
of an electromagnetic 3D-Sinai billiard was analyzed ac¬ 
cording to standard methods of Random Matrix and Pe¬ 
riodic Orbit Theory. Spectral correlations are shown to 
be totally consistent with the predictions of the GOE af¬ 
ter the systematic extraction of the family of 3D bbos 
which dominate not only the length spectrum of the bil¬ 
liard but also lead to dramatic deviations from Gaus¬ 
sian characteristics in the measures for long-range corre¬ 
lations, E 2 and A 3 . 
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